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We revisit on rational solution of Sasa-Satsuma equation, which can be used to describe evolution 
of optical field in a nonlinear fiber with some high-order effects. We find a striking dynamical 
process which involves both modulational instability and modulational stability regimes, in contrast 
to the rogue waves and W-shaped soliton reported before which involves modulational instability 
and stability respectively. It is demonstrated that stable W-shaped solitons can be generated from 
a weak modulation signal on continuous wave background. This provides a possible way to obtain 
stable high-intensity pulse from low-intensity continuous wave background. 
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I. INTRODUCTION 

Recently, rational solutions of nonlinear Schrddinger 
equation (NLSE) have been paid much attention, since 
rational solution of NLSE can be used to describe rogue 
wave (RW) dynamics in many different physical sys¬ 
tems [ll-Q- Among these different systems, optical 
fiber play an important role in experimental observa¬ 
tions for its well-developed intensity and phase modu¬ 
lation techniques. The experimental studies in nonlin¬ 
ear fiber showed that the simplified NLSE can well de¬ 
scribe the dynamics of localized waves, which only con¬ 
tains the group velocity dispersion (GVD) and its coun¬ 
terpart, namely, self-phase modulation (SPM). But for 
ultrashort pulses whose duration is shorter than lOOfs, 
which is tempting and desirable to improve the capac¬ 
ity of high-bit-rate transmission systems, the nonlinear 
susceptibility will produce higher-order nonlinear effects 
like the Kerr dispersion (i.e., self-steepening) and the de¬ 
layed nonlinear response except for SPM, and even the 
third-order dispersion (TOD). These are the most general 
terms that have to be taken into account when extend¬ 
ing the applicability of the NLSE d, 0 . With these ef¬ 
fects, the corresponding integrable equation was derived 
as Sasa-Satsuma(S-S) equation M- 

The linearized stability analysis for the S-S equation 
in Ref. m suggested that there are both modulational 
instability(MI) and modulational stability(MS) regimes 
for low perturbation frequencies on the continuous wave 
background(CWB). In the MI regime, the rational solu¬ 
tions are obtained in [I^, Il3| : which correspond to RW 
excitation. It is demonstrated that the high-order effects 
could make the RW twisted and the rational solution of 
the S-S equation had distinctive properties in contrast 
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to that of the well-known NLSE. In the MS regime, we 
presented rational solutions in the MS regime, which cor¬ 
respond to W-shaped soliton excitation [ij. It is shown 
that the profile of W-shaped soliton depends on the back¬ 
ground frequency. The results suggested that not all ra¬ 
tional solution of nonlinear partial equation correspond 
to RW dynamics. It should be noted that there are two 
critical points for background frequency on boundary of 
MI and MS regime. We expect that there should be some 
new dynamics on the critical points. 

In this paper, we revisit on rational solution of S-S 
model on the critical background frequencies. We find 
two W-shaped solitons can be generated from a weak 
modulation signal on continuous wave background. The 
striking process should involves both modulational insta¬ 
bility and modulational stability regimes, in contrast to 
the rogue waves and W-shaped soliton reported before 
which involves modulational instability and stability re¬ 
spectively. 

II. TWO W-SHAPED SOLITONS GENERATED 
FROM A WEAK MODULATION ON CWB 

For ultrashort pulses whose duration is shorter than 
lOOfs, the optical field [i.e., E{t,z)] in a nonlinear fiber 
with some high-order effects can be described by the S-S 
model [T^ : 

iE^ + ^Ett + \E\^E 

+ie[Ettt + 3(|EpE)t + 3|EpEt] = 0, (1) 

where e is an arbitrary real parameter to scale the in¬ 
tegrable perturbations of the NLSE and it can change 
the strength of the high-order effects conveniently. The 
last three high-order terms correspond to the TOD, self- 
steepening, and delayed nonlinear response respectively 
[^ flR - flTj l. Raman effect TnElEl"^ {Tr is a real constant) 
is not considered in the model, since the term can be 
much smaller than Kerr dispersion effect, and it makes 
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FIG. 1: (Color online) The evolntion of IV'I^ [given by Eq. 
®1 corresponding to the process that one weak modnlation 
signal evolves into two W-shaped solitons. The parameters 
are a = 71 = 72 = 1 . 
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FIG. 2: (Color online) The corresponding spectrnm evolution 
of the process illustrated in Fig. 1. It is seen that there is an 
interference pattern on spectral distribution with ijj' = uj — ojo- 


the pulse energy non-conservation (this usually makes 
the model non-integrable). When e = 0, the S-S equa¬ 
tion reduces to the standard NLSE which has only the 
terms describing the lowest order dispersion and self¬ 
phase modulation. The investigations on the S-S equa¬ 
tion indicated that the nonlinear waves in nonlinear fiber 
with the high-order effects are much more diverse than 
the ones in simplified NLSE (Tbl - f^ . Here, we revisit on 
rational solution of the S-S model. 

Eor convenience, let tpiT, Z) = 

E{t,z)exp[--^{t- [U, where T = t - and 
Z = z, then a simplified equation equivalent with Eq. 
0 is obtained as i/fz -I-e(V'TTT + = 0. 

The evolution direction of nonlinear waves in this 
equation is slant, which make us inconvenient to observe 
the dynamical behaviors. To overcome this difficulty, 
we introduce the coordinates transformation, namely, 
X = a (T — and rj = —^^a^eZ, then the 

simplified equation can be written as 


wave near r] = —1.5. In particular, this large localized 
wave splits into two quasi-stabilized nonlinear waves after 
r] = —1.5. The profiles of them can be observed clearly 
with the help of making cut plot at a certain location ry, 
which show a nice W-shaped soliton configuration simi¬ 
lar to the one reported in Ref. [l^. It is found that the 
higher peak of the two waves decreases with increasing 
the propagation distance rj while the lower one increases. 
However, after a certain distance, the changing rates of 
them will tend to zero, namely, the nonlinear localized 
waves tend to solitons whose shapes are kept well. This 
process could suggest one way to excite quasi-stabilized 
W-shaped soliton through adding small modulation on 
the continuous background, which has significant and po¬ 
tential values in the field of soliton application. 

The above dynamical process can be faithfully de¬ 
scribed by the following analytical solution on the CWB 
with frequency ojq = ±i, which can be obtained by 
adopting Darboux transformation method. 


11^3 
= - 
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The coordinates transformation only lead to a trivial 
phase multiplication from Eq. ([2]) to Eq. (ED and thus 
the identical dynamics are kept well. Subsequently, we 
will discuss the nonlinear dynamics of Eq. m instead 
of Eq. ([T]) without losing generality. On the CWB 
'0o(x,??) = a exp [iojoX - + ^u}^)r]] with arbi¬ 

trary ujQ, it is convenient to study various localized waves 
such as RW with two peaks [H, [l^ and rational W- 
shaped soliton M- Surprisingly, in the case of wq = ± 5 , 
we find a dynamics process that one weak modulation 
signal on the background can evolve into two stable non¬ 
linear waves as shown in Fig. I. This process is very 
different from the dynamics of nonlinear waves for the 
same model reported in Refs. [iMl. From Fig. 1 we 
see that the propagating signal at the locations before 
rj = —2 keeps small while it evolves to a large localized 
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with 


H = iLi(x,r/) + (2zy3-12zx-12l72-6)7i 

-l-[I 8 zx^ - (I8zx/3 -I- 18)x - 12 z ?7 -I- 6x/3 -I- 6 z ]72 
+ (12fx-8zy3-12)7|. 


M = Mi(x, Z7) +127?+ (12x'-16V3x +20)72" 

+ [(16^3 - 24x) 72 + 24?? - 12x" + lOv^x - 8]7i 
-b[12x^ - 2AV3x^ + (48 - 24?7)x + IdVdv - 8\/3]72, 

where 71 and 72 are two real constants. 7Li(x,»?) = 
6 zx^ — (9?-\/3-|-9)x" + {dVi + di— 12 z 77 )x+ (2z\/3 — Q)ri + 
2i\fd — 6 and Mi(x, ??) = 3x^ — 8\/3x^ + (24 — 12??)x" + 
(16-\/3?? — 8-\/3)x + 4 — 8 ??-|- 12??^. When a = 71 = 72 = 1, 
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FIG. 3: (Color online) (a) The MI gain distribution on back¬ 
ground frequency and perturbation frequency plane. The 
crossing points of black lines and white line are the criti¬ 
cal points between MI and MS for resonant perturbation fre¬ 
quency. (b) The cut plot of MI gain on uj' = 0. It is shown 
clearly that the critical points are at ±1/2. 


the dynamics of solution ([3|) corresponds to the one in 
Fig. 1. For a ^ \ and other values of 71 ^ 2 , this solu¬ 
tion demonstrates the identical dynamical process except 
some trivial shifts on the temporal-spatial plane. Per¬ 
forming asymptotic analysis on the exact solution, we 
find that the two W-shaped solitons will have identical 
profiles when 77 —^ 00 . The hump intensities and valley 
intensities of them tend to be 4 and | times the CWB 

O 

intensity value, respectively. This analysis indicates that 
the two W-shaped localized waves indeed evolve into soli¬ 
tons. 

As shown in Ref. [l^, the spectrum of W-shaped 
soliton corresponds to one stable broad spectrum pulse. 
This means that the two W-shaped solitons correspond 


to two broad spectrum pulses, which can be proved by 
making Fourier transformation on two asymptotic soli¬ 
ton forms of the exact solution. Since there are two 
W-shaped solitons generated from a small modulation 
signal, one can expect emerging two stable triangular 
spectrum and interference fringes. The spectrum evo¬ 
lution of the dynamical process illustrated in Fig. 2 
is obtained exactly through the Fourier transformation 

The back¬ 
ground is infinity and it is integrated to be S(uj — cuq) 
(here the background frequency ojq = 5 ), therefore we 
can eliminate the (5-function by taking the transforma¬ 
tion on the varying part. These properties have been 
checked (see Fig. 2) and the expected interference pat¬ 
tern on the spectral distribution has been observed. 

Now we try to understand the above process for which 
two W-shaped solitons excited from a small perturbation. 
To this end, we first perform the linear stability analysis 
(LSA) on the CWB v) by adding the Fourier modes 
with small amplitudes [lH, i.e., /+exp [iw'/x — iFr;)] ± 
/I exp [—ioj'{x ~ K*ri)\. The MI gain spectrum on back¬ 
ground frequency and perturbation frequency space is 
shown in Fig. 3(a). The results show that there are 
two different regimes for perturbations on the continuous 
background, namely, MI and MS regimes. The critical 
points for background frequency are Wq = as shown 
in Fig. 3(b). For |a;o| < there exists a MS regime 
for which the perturbations with arbitrary frequencies 
are all stable on the CWB. In this regime, stable ratio¬ 
nal W-shaped soliton solution has been obtained in Ref. 
M- For I Wo I > there exists a nonzero MI regime 
in which the perturbations with frequencies belonging 

to [wq — ^4 — uJq^, Wo ± ^4 — W(/^], and the maximum 
gain of MI will emerge on the “resonant” line for which 
the perturbation frequency is equal to the CWB’s. The 
resonant line in MI regime usually corresponds to RW 
excitation. In this regime, RW solution which possesses 
double humps has been given in Refs. [H, [I3|. Notably, 
the two W-shaped solitons here growing from a small 
perturbation should involve both two regimes. This pro¬ 
cess is found on the CWB with the critical frequencies 
Wo = ±^ which is located on the boundary between the 
MS and MI regime. We can qualitatively know that the 
first stage of the process is in the MI regime for which the 
small modulation signal grows to one high hump and two 
deep valleys and then it splits. The second stage of the 
process should evolve in the MS regime, which makes the 
split localized waves stable and evolve into two W-shaped 
solitons. This indicates that MI and MS can coexist on 
the critical background frequency given by the linear sta¬ 
bility analyze. The MI analysis on Eq. © shows that 
the MS regime for low-frequency perturbations does not 
exist when neglecting the high-order effects. This implies 
that these high-order effects play important roles in the 
dynamical process. 

To see the feasibility of process for W-shaped soli¬ 
tons generation, we numerically evolve Eq. m from 
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FIG. 4: (Color online) (a) The numerical simulation for the 
dynamics of |'0(Xi’7)l^ evolved from an initial given by the 
exact solution at r; = —3. (b) The corresponding evolution 
given by the exact analytical result. It is shown that they 
agree well with each other, which indicates that the exact 
signal evolution can survive even with numerical deviations. 
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FIG. 5: (Color online) The profiles of the nonlinear waves 
at different distances. The solid gray curves denote the an¬ 
alytical solutions while the dashed ones correspond to the 
numerical results. 


an exact initial signal with ry = —3 by adopting the 
split-step Fourier method and the results for are 
shown in Fig. 4(a) compared with the analytical so¬ 
lution [see Fig. 4(b)]. Figure 5 demonstrates the pro¬ 
files comparison. It is shown that the simulation results 
agree well with the analytical ones even with numerical 
deviations. A bit differences for the humps’ locations 
and the intensities of the humps is caused by the peri¬ 
odic boundary condition induced reflection effects. The 
numerical simulation from the initial small signal with 
some additional noise still evolves to two W-shaped soli- 
tons. This results suggest that the exact process might 
be realized experimentally. It must be mentioned that 
the W-shaped solitons discussed here are grown from a 
small modulation signal on the CWB, which are differ¬ 
ent from the RW and W-shaped solitons reported in ear- 
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FIG. 6: (Color online) The deviations between |'!/'i|^ and 
|'!/) 2 p. The numerical evolutions are from the initial condi¬ 
tions that V’i(x, -3) = V'(X. -3) and 'i/’ 2 (x, -3) = i/’(x, -3) + 
0.01 Random[x\. It is shown that the deviations are lower 
than 0.04. 

her studies [l^ - [T3| . Furthermore, we test the stability 
of the generated W-shaped solitons with adding a white 
noise 0.01 Random(x) ix ^ [“I;!])- an example, we 
show the deviations between the numerical result with 
and without the noise in Fig. 6. We can see that the 
solitons are stable against the noise even with the Ra¬ 
man term. 


III. CONCLUSION AND DISCUSSION 

In summary, we present a striking dynamical process 
in which two W-shaped solitons can be generated from 
a weak modulation signal on the CWB. The numerical 
simulations indicate that the dynamical process is robust 
against small noise or perturbations. The underlying 
mechanisms of the process has been discussed qualita¬ 
tively based on the MI analysis. The process involves MI 
to MS autonomously. This indicates that MI and MS 
can coexist on the critical background frequency given 
by LSA. This is a novel nonlinear dynamical behavior, 
which could exists in many other nonlinear systems, and 
even nonlinear coupled systems. Notably, the critical fre¬ 
quency located on the boundary between the MS and MI 
regimes is obtained by LSA. The analysis predicts there 
is no MI on the critical frequency, but we demonstrate 
that both MI and MS play key roles in the dynamical 
process. This means that LSA could not stand precisely 
on the critical line. The explicit underlying mechanism 
for the dynamical process still needs new analysis tech¬ 
niques. 

Actually, many recent experiments for nonlinear local¬ 
ized waves including Kuanetsov-Ma breather, Akhmediev 
breather, and Peregrine RW in nonlinear fibers 
have provided a good platform to investigate complicated 
dynamics of nonlinear localized waves on CWB. These 
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experimental works suggested that the exact analytical 
solutions for the simplified NLSE can well describe the 
evolution of optical fields even with non-ideal initial con¬ 
ditions. However, the ideal initial conditions is impor¬ 
tant to study the NLSE with high-order effects and our 
present work could provide the initial conditions for gen¬ 
erating W-shaped solitons in experiment from weak mod¬ 
ulation on CWB. This has been checked numerically even 
for the initial conditions with small noises. In particular, 
our theoretical results can be used to qualitatively un¬ 
derstand the experiment in a nonlinear fiber with TOD 
and self-steepening effects [2^. In this experiment one 
Peregrine RW split into two lower-amplitude localized 
solitons and the pulse splitting process was explained by 


high-order MI based on a high-order Akhmediev breather 
solution of the simplified NLSE [ 2 ^. Here, we provide 
another perspective for understanding the pulse splitting 
process and predict the quasi-stability of the splitting 
pulses with some certain high-order effects. 
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